arXiv:1507.07353v2 [astro-ph.GA] 18 May 2016 


Astronomy & Astrophysics manuscript no. tf_accepted ©ESO 2016 

May 19, 2016 


Disk galaxy scaling relations at intermediate redshifts 
I. The Tully-Fisher and velocity-size relations * 

Asmus Bohm 1 and Bodo L. Ziegler 2 


1 Institute for Astro- and Particle Physics, Technikerstrasse 25/8, 6020 Innsbruck, Austria, e-mail: asmus. boehm@uibk. ac. at 

2 Institute of Astronomy, Tiirkenschanzstrasse 17, 1180 Vienna, Austria. 

Received ; accepted 


ABSTRACT 

Aims. Galaxy scaling relations such as the Tully-Fisher relation (between maximum rotation velocity V max and luminosity) and the 
velocity-size relation (between V max and disk scale length) are powerful tools to quantify the evolution of disk galaxies with cosmic 
time. 

Methods. We took spatially resolved slit spectra of 261 field disk galaxies at redshifts up to z ~ 1 using the FORS instruments of the 
ESO Very Large Telescope. The targets were selected from the FORS Deep Field and William Herschel Deep Field. Our spectroscopy 
was complemented with HST/ACS imaging in the F814W filter. We analyzed the ionized gas kinematics by extracting rotation curves 
from the 2-D spectra. Taking into account all geometrical, observational and instrumental effects, these rotation curves were used to 
derive the intrinsic V max . 

Results. Neglecting galaxies with disturbed kinematics or insufficient spatial rotation curve extent, V max could be robustly determined 
for 124 galaxies covering redshifts 0.05 < z < 0.97. This is one of the largest kinematic samples of distant disk galaxies to date. 
We compared this data set to the local 6-band Tully-Fisher relation and the local velocity-size relation. The scatter in both scaling 
relations is a factor of ~ 2 larger at z » 0.5 than at z ~ 0. The deviations of individual distant galaxies from the local Tully-Fisher 
relation are systematic in the sense that the galaxies are increasingly overluminous towards higher redshifts, corresponding to an over¬ 
luminosity A M b = —(1.2 ± 0.5) mag at z = 1. This luminosity evolution at given V max is probably driven by younger stellar populations 
of distant galaxies with respect to their local counterparts, potentially combined with modest changes in dark matter mass fractions. 
The analysis of the velocity-size relation reveals that disk galaxies of a given V max have grown in size by a factor of ~ 1.5 over the past 
~ 8 Gyr, likely via accretion of cold gas and/or small satellites. Scrutinizing the combined evolution in luminosity and size, we find 
that the galaxies which show the strongest evolution towards smaller sizes at z ~ 1 are not those which feature the strongest evolution 
in luminosity, and vice versa. 
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1. Introduction 


Observational studies of galaxy evolution have made great 
progress in the past decade. In the course of some projects, red¬ 
shifts and spectral energy distributions of several 10 4 galaxies 
at significant cosmological look-back times were gained — two 
examples are the VVDS (Le Fevre et al. 120051) and zCOSMOS 
(Lilly et al. 120071 ) surveys. Other projects focused on smaller 
samples to conduct detailed studies of scaling relations which 
link fundamental structural and kinematical parameters of galax¬ 
ies. In the case of spirals, the most famous scaling relation is 
the Tully-Fisher Relation (TFR, Tully & Fisher 11977 ) which 
relates the luminosity to the maximum rotation velocity V max . 
Equivalents of the classical, i.e. optical TFR were later estab¬ 
lished, showing that also the stellar mass M, (e.g. Bell & de 
Jong 1 2001 ) or total baryonic mass (i.e. stars and gas, e.g. Mc- 
Gaugh et al. 120001 correlate with V max . In that sense, the optical 
TFR is only a variant of a more fundamental relation between the 
baryonic and the dark matter content of disk galaxies. Since dark 
matter dominates the total mass budget, V max can be used to es¬ 
timate the dark matter halo mass (e.g. Mo, Mao & WhiteJT998 ). 


* Based on observations with the European Southern Observatory 
Very Large Telescope (ESO-VLT), observing run IDs 65.0-0049, 66. A- 
0547, 68^-0013, 69.B-0278B, 70.B-0251A and 081.B-0107A. 


Besides luminosity, stellar or baryonic mass, also the disk scale 
length is correlated with V max (e.g. Mao, Mo & White 119981 ). 
this is referred to as the rotation velocity - size relation (VSR). 
The three parameters V max , size and luminosity span a parame¬ 
ter space in which disks populate a two-dimensional plane with 
small intrinsic scatter. E.g., Burstein et al. d 19971 ) utilized this pa¬ 
rameter space characterizing size via the effective radius; Koda 
et al. (120001 ) carried out a similar investigation using the disk 
isophotal radius at p/ = 23.5 mag. The distribution of late-type 
galaxies within a plane in this parameter space is similar to the 
fundamental plane valid for early-type galaxies (e.g. Dressier 

et al. [19571 . 

The first observational attempts to construct the optical 
TFR of distant spirals and to quantify their evolution in lu¬ 
minosity were made almost two decades ago, e.g. by Vogt 
et al. d 19961 ) and Rix et al. (1 1 9971 . For many of the following 
years, there have been discrepant results from different studies 
on a possible evolution of the Tully-Fisher relation with cos¬ 
mic time in zero point or slope. Regarding the //-band TFR, 
Vogt (120011 did not find any evolution up to redshifts z ~ 1, 
while e.g. Bohm et al. (120041 . Bamford et al. (120061 or Fer¬ 
nandez Lorenzo (120101 found that disk galaxies at z ~ 1 were 
brighter by A Mg ~ -1 mag for a given maximum rotation ve¬ 
locity. Bohm et al. also discussed a possible slope change with 
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cosmic time, in the sense that the luminosity evolution of low- 
mass spirals was stronger than that of high-mass ones. Weiner 
at al. (l2006t >. on the other hand, found the opposite evolution, 
i.e. a stronger brightening in high-mass disk galaxies. Bohm & 
Ziegler ( 120071 1 showed that a strong evolution of the TFR scatter 
could mimic an evolution in TFR slope due to selection effects. 
Kassin et al. (120071 1 established that even galaxies with kinematic 
disturbances — which usually do not follow the classical TFR — 
obey a remarkably tight correlation when non-ordered motions 
are taken into account. To this end, these authors introduced the 
parameter S 0.5 which combines V max and gas velocity dispersion, 
finding a constant slope over the epoch 0.1 < z < 1.2. 

Only more recently, there has been growing consent that the 
local TFR slope holds at least up to redshifts about unity, for 
all variants of the TFR. Some of these studies were — like all 
of the projects mentioned above — based on slit spectroscopy 
(e.g. Fernandez Lorenzo et al. 12010 1 Miller et al. 12011 ). others 
made use of Integral Field Units (IFUs, e.g. Flores et al. 12006 1 
Puech et al. 12008 1). Due to the time-expensive approach, samples 
constructed using IFUs usually are smaller than slit-based ones, 
the gain is a direct observability of (at least part of) the two- 
dimensional rotation velocity field. Mismatches between photo¬ 
metric and kinematic center, or photometric and kinematic posi¬ 
tion angle, can only be detected with this kind of data. 

Towards high redshifts z > 1, the usage of IFUs is a ne¬ 
cessity, since interaction and merger events were much more 
frequent at these epochs. In effect, large fractions of high-z 
star-forming galaxies feature complex or disturbed kinematics. 
Several surveys have been conducted with the adaptive optics- 
assisted SINFONI instrument of the ESO Very Large Telescope: 
AMAZE (Maiolino et al. 120081) . SINS (Forster-Schreiber et 
al. 120091 ) and MASSIV (Epinat et al. 120091) . Based on data from 
MASSIV, e.g. Vergani et al. (120121) found only a small increase 
in stellar mass (~ 0.15 dex on average, depending on the refer¬ 
ence sample at z = 0) between z ~ 1.2 and the present-day uni¬ 
verse, at fixed V max . Gnerucci et al. (1201 il l investigated 11 disks 
at z ~ 3 from the AMAZE data set, finding stellar masses smaller 
than locally by ~ 1 dex. However, even for this sample of regu¬ 
larly rotating disks, a very large scatter was observed, and the 
authors concluded that the TFR is not yet established at that cos¬ 
mic epoch. Cresci et al. ( 120091 ). on the other hand, found a much 
smaller evolution of 0.41 dex in log M, since z ~ 2.2 using SINS 
data. 

A potential environmental dependence of the TFR has been 
subject to many studies. Several authors found that field and 
cluster samples have the same TFR slope, but the scatter is in¬ 
creased in dense environments (e.g. Moran et al. 120071 Bosch et 
al. [2013b ). This probably is induced by cluster-specific interac¬ 
tion processes. Tidal interactions between close galaxies can in¬ 
crease the star formation rate (e.g. Lambas et al. 12003 ). whereas 
interactions between the interstellar medium and the hot intra¬ 
cluster medium (ram-pressure stripping), can push gas out of a 
disk galaxy and in the extreme case totally quench star forma¬ 
tion (e.g. Quilis et al. 2000). These mechanisms lead to a larger 
range in luminosities at given V max and increase the fraction of 
perturbed gas kinematics (e.g. Bosch et al. 12013a] ). In turn, V max 
measurements in dense environments carry larger systematic er¬ 
rors. All these effects are likely to contribute to the larger TFR 
scatter found in clusters. However, the situation changes when 
field and cluster samples are matched in rotation curve quality, 
i.e. when galaxies are rejected that have perturbed kinematics 
due to cluster-specific interactions. It was found that the distribu¬ 
tions of field and cluster samples then are very similar in Tully- 


Fisher space (e.g., Ziegler et al. , 20031 Nakamura et al. 12006; , 
Jaffeet al. l20TT1) . 

Concerning numerical simulations, spiral galaxies often had 
too low an angular momentum compared to observed ones (e.g. 
Steinmetz & Navarro[19991). Only more recently, it became fea¬ 
sible to simulate galaxies which over a broad mass range agree 
with the observed, local TFR, by including recipes for internal 
physics such as supernova feedback as well as external processes 
such as the ultraviolet background (e.g. Governato et al. 120071 ). 
Dutton et al. (1201 11 1 used combined /V-body simulations and 
semi-analytic models to predict the evolution of several scal¬ 
ing relations up to redshift z = 4. In these simulations, disks 
at z = 1 are, at given Vm ax , brighter by -0.9 mag in the /I-band 
and smaller by ~ 0.2 dex than their local counterparts. Based on 
cosmological N - b o d y / h y d ro d y n am i c a I simulations, also Porti- 
nari & Sommer-Larsen (120071) found a fi-band brightening by 
-0.85 mag at z - 1. 

Observational studies of the evolution of the VSR are rela¬ 
tively scarce. Puech et al. (120071) found no change in disk sizes 
at given V max between z ~ 0.6 and z = 0. Vergani et al. (120121) 
reported only a small increase of 0.12 dex in half-light radius 
since z ~ 1.2. Towards higher redshifts, a stronger evolution was 
found: the sample presented by Forster-Schreiber et al. (120091) 
yields an increase in size by a factor of ~ 2 between z ~ 2 and 
locally (see Dutton et al. 1201 II . However, this value might be an 
underestimate due to the computation of disk sizes based on Her 
half-light radii, as shown by Dutton et al. 

In this paper, we will use the Tully-Fisher and velocity-size 
relations to investigate the evolution of disk galaxies in luminos¬ 
ity and size since redshifts z ~ L Note that the main driver of 
our project is not a complete census of the disk galaxy popula¬ 
tion during these cosmic epochs, but a detailed look at only the 
virialized and undisturbed disks. Only this allows to use scaling 
relations like the Tully-Fisher without the impact of kinematic 
biases. The paper is organized as follows: In Sect. [2] we outline 
the selection and observation of our sample. Sect. [3] briefly de¬ 
scribes the data reduction, in Sect. [4] we construct and analyze 
the intermediate-redshift scaling relations. Sect, [^comprises the 
discussion and Sect.[6]summarizes our main results. 

In the following, we assume a flat concordance cosmology 
with Da = 0.7, D,„ = 0.3 and Ho = 70kms _1 Mpc -1 . All mag¬ 
nitudes are given in the Vega system. 

2. Sample selection and observations 

For the selection of our spectroscopic targets, we relied on two 
multi-band photometric surveys: the FORS Deep Field (FDF; 
Heidt et al. 120031) and the William Herschel Deep Field (WHDF; 
Metcalfe et al. 12001 1). These comprise deep imaging in the filters 
U , B, g , R, /, J , K (FDF) and U , B, R, /, H , K (WHDF). The filter 
set used in the FDF photometry is very similar to the Johnson- 
Cousins system, while the WHDF photometry is based on Harris 
filters; we transformed these magnitudes to the Johnson-Cousins 
system via synthetic photometry. 

We applied the following criteria to construct the sample: 

1. total apparent brightness R < 23 mag; 

2. star-forming spectral energy distribution, based on a pho¬ 
tometric redshift catalog for the FDF targets (Bender et 
al. 12001) or color-color diagrams for galaxies in the WHDF, 
for which no photometric redshifts were available. For the 
latter candidates, we adopted the evolutionary tracks pre¬ 
sented by Metcalfe et al. (12001) : 
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3. disk inclination angle i > 30° to avoid face-on disks and 

ensure sufficient rotation along line-of-sight; 

4. misalignment angle 6 < 15° between apparent major axis 

and slit direction to limit geometric distortions of the ob¬ 
served rotation curves. 

Note that no selection on morphological type nor emission line 
strength was used. 

The spectroscopic data were taken between September 2000 
and October 2008 with the FORS 1 & 2 instruments of the 
VLT. In total, 261 disk galaxies were observed. All runs except 
the one in 2008 were carried out in multi-object spectroscopy 
(MOS) mode with straight slitlets perpendicular to the direction 
of dispersion. The observations in 2008 made use of the Mask 
Exchangeable Unit (MXU) with tilted slits, allowing to accu¬ 
rately place them along the apparent major axes and achieve a 
misalignment angle 6 « 0°. Slit tilt angles 6 were limited to 
|0| < 45° to ensure a robust sky subtraction and wavelength cal¬ 
ibration. We used a fixed slit width of 1.0 arcsec which resulted 
in a spectral resolution of R « 1200 (grism 600R) and a spa¬ 
tial scale of 0.2 arcsec/pixel for observations taken before 2002 
when the FORS CCD was upgraded. The upgrade led to an in¬ 
creased sensitivity at wavelengths A > 7000 A, a lower readout 
noise and much higher readout speed. The data taken after this 
upgrade feature R ~ 1000 (grism 600RI) and a scale of 0.25 
arcsec/pixel. The total integration time for all MOS setups was 
2.5 h. Only for the MXU observations in 2008, we used a much 
longer integration time of 10 h total per target. This latter run was 
designed to particularly extend our sample at low luminosities, 
and included galaxies down to an apparent A 1 -band brightness of 
R ~ 24. Except for fill-up targets, all galaxies from this run have 
a /(-band absolute magnitude Mb > -19. Across all our spectro¬ 
scopic campaigns, seeing conditions ranged from 0.42 arcsec to 
1.20 arcsec FWHM, with a median of 0.76 arcsec. 

To determine structural parameters such as disk inclination 
and scale length, bulge-to-disk ratios etc., we took HST/ACS 
images with the F814W filter (similar to the Cousins /-band). 
The 6x6 arcmin 2 sky areas of the FDF and WHDF were cov¬ 
ered with a 2 X 2 mosaic each, using the Wide Field Camera 
(0.05 arcsec/pixel) with one orbit per pointing and total integra¬ 
tion times of 2360 s (FDF) and 2250 s (WHDF), respectively. 

3. Data reduction 

All reduction steps were carried out on the extracted 2-D spectra 
of each exposure and each target individually. The reduction in¬ 
cluded bias subtraction, cosmics removal, flatfielding, correction 
of spatial distortions, wavelength calibration and sky subtraction. 
Only the final, fully-reduced 2-D spectra were co-added with a 
weighting factor according to the seeing conditions during spec¬ 
troscopy. For wavelength calibration, the dispersion relation was 
fitted with a third-order polynomial row by row; the median rms 
of these fits was 0.04 A. 

For the data taken with tilted slitlets in 2008, we took a differ¬ 
ent approach and applied the improved sky subtraction method 
by Kelson (120031) . It avoids the problems arising from the trans¬ 
formation of a regular pixel grid (raw spectrum) to an irreg¬ 
ular pixel grid (rectified, wavelength-calibrated spectrum). To 
this end, the night sky emission is fitted and removed after bias 
subtraction and flatfielding, but before the distortion correction 
and wavelength calibration. Note that the difference in final data 
quality between the “classical" method, where sky subtraction is 
carried out as the last of all reduction steps, and the one follow¬ 
ing Kelson is only marginal when slitlets are oriented perpendic- 



Fig. 1. Redshift distribution of all galaxies in our survey (solid line, 
238 galaxies; one fill-up target at z = 1.49 is omitted in this plot) and 
our kinematic sample, i.e. all galaxies used in our kinematic analysis 
(hashed histogram, 124 galaxies). The selection criterion for the latter 
sub-sample is detailed in Sect. 14. 31 

ular to the direction of dispersion (as is the case for our FORS 
data taken in MOS mode). However, for data gained with tilted 
slits (MXU mode), night sky line residuals are strongly reduced 
with the Kelson approach (Kelson 2003! ). 

Regarding the ACS imaging, the standard pipeline was used 
to carry out bias subtraction, flatfielding and distortion correc¬ 
tion. We applied a filtering algorithm to finally combine the ex¬ 
posures of each pointing and remove the cosmics. 

4. Analysis of the scaling relations 

Out of the 261 disk galaxies for which we obtained spectra, red- 
shifts could be determined for 238 galaxies. These objects range 
from z = 0.03 to z = 0.97 (omitting an outlier at z = 1.49 which 
was a fill-up target and did not yield a V max value) with a median 
of (z) = 0.43. We show the redshift distribution in Fig.Q] The gap 
around z ~ 0.5 can be attributed to the galaxies stemming from 
the FDF. This gap is not a “redshift desert" due to constraints in 
our spectroscopic setup; it is also present in the distribution of 
the FDF photometric redshifts, while it is absent in the redshift 
distribution of the WHDF galaxies only. Most probably it thus is 
physical and a result of cosmic variance — both the FDF and the 
WHDF are deep surveys with a relatively small field-of-view of 
» 40 arcmin 2 each. 

4 . 1. Absolute magnitudes 

We opt to construct the /(-band Tully-Fisher relation in this study 
to be sensitive for recent or ongoing star formation. A big advan¬ 
tage for this is the multi-band imaging at hand closely match¬ 
ing the rest-frame /(-band, for any galaxy redshift in our sam¬ 
ple. To determine the /(-band luminosities Mb, we used the ob¬ 
served filter X obs best matching the rest-frame /(-band to de¬ 
rive the /-correction via synthetic photometry and compute the 
transformation X obs —» B res t- In the FDF, we used the trans¬ 
formation B obs —> Brest at redshifts z < 0.25, g 0 bs —» Brest at 
0.25 < z < 0.55, R 0 bs —> B r est at 0.55 < z < 0.85 and / obs —> B res t 
at z > 0.85. For the galaxies stemming from the WHDF, we uti¬ 
lized B 0 bs —» Brest at z < 0.3, B obs —» B res t at 0.3 < z < 0.7 and 
fobs —» B re st at z > 0.7. Thanks to this approach, the /-correction 
is only weakly depending on the spectral energy distribution of 
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a given galaxy: if the spectral classification would be wrong by 
AT = 2 (which corresponds to the spectrum of, for example, an 
Sa galaxy mistakenly classified as type Sb), the resulting system¬ 
atic k-correction error would be cr> <0.1 mag across the whole 
redshift range covered by our data. 

We corrected for intrinsic absorption A# due to the dust disk 
following the approach by Tully et al. (I1998i >. This formalism is 
inclination- and V max -dependent: disks which are observed more 
edge-on have a higher extinction than more face-on ones, and 
more massive disks, i.e. galaxies with higher L max , have a higher 
extinction than less massive ones. 

To summarize, the absolute /(-band magnitudes were com¬ 
puted as 

M b = m x — A s x — DM — k B — A B , (1) 

where nix is the total apparent brightness in filter X , A s x is the 
Galactic absorption for filter X, corrected using the maps by 
Schlegel et al. ( il9981 >. DM is the distance modulus, k H is the k- 
correction and A b is the correction for intrinsic dust absorption 
in rest-frame B. 

4.2. Structural parameters 

We derived the structural parameters of the galaxies — disk in¬ 
clination, disk scale length rj, bulge-to-total ratio B / T etc. — 
on the HST/ACS images using the GALFIT package by Peng 
et al. (120021 1. It allows to fit multiple two-dimensional surface 
brightness profiles simultaneously to the galaxy under scrutiny 
as well as to neighboring objects (the importance of such simul¬ 
taneous fits is discussed, e.g., in HauBler et al. 120071 . GALFIT 
requires an input Point Spread Function (PSF) for the convo¬ 
lution of the model profiles. We constructed PSFs for the FDF 
and WHDF separately from ~ 20 unsaturated stars in each field. 
Both have a FWHM of 0.12arcsec, corresponding to a spatial 
resolution of ~ 0.7 kpc at z = 0.5. 

The surface brightness profiles of all galaxies in our sample 
were fitted with two different setups: i) a single Sersic profile 
with free index « ser or ii) a two-component model with an expo¬ 
nential profile for the disk and a Sersic profile with fixed index 
/z s er = 4 for the bulge. The best-fit parameters from method i) 
were used as initial guess values for method ii). 

All fit residuals were visually inspected, and in a few cases, 
constraints on parameters like e.g. bulge effective radius were 
necessary to avoid a local x 2 minimum in the fitting process. 
For our analysis, we mostly used the disk parameters from the 
bulge/disk decomposition, i.e. method ii). Only in a few cases 
where the two-component fit of an evidently bulgeless disk was 
not converging, we kept the parameters from the single Sersic fit. 
We stress that for the analysis presented here, the most impor¬ 
tant parameters are the inclination i, position angle 9 and scale 
length /■(] of the disk. These showed only small differences be¬ 
tween the two fitting methods. This is mainly because the vast 
majority of the FDF/WHDF disks have only small bulges or 
even no detectable bulge at all; the median bulge-to-total ratio 
is {B/T) = 0.06. 

GALFIT only returns random errors on the best-fit param¬ 
eters. These are very small (<1 %) throughout our sample. To 
gain a more realistic estimate of the systematic errors on r&, we 
can rely on our own previous analysis of HST/ACS images us¬ 
ing GALFIT in Bohm et al. (120131) . In that work, we investigated 
the impact of an Active Galactic Nucleus on the morphologies of 
host galaxies at redshifts 0.5 < z < 1.1 as quantified with GAL¬ 
FIT. For a negligible central point source, we found a typical 


systematic error of 20 % on galaxy sizes (see Fig. 7d in Bohm 
et al. [20131 . This value hence represents the systematic size er¬ 
ror for galaxies with the light profiles of pure disks or disks with 
only weak bulges; this is the case for the vast majority of galaxies 
in our sample. We therefore adopt this error on r d in the follow¬ 
ing. 

It is well-known that the observed disk scale length depends 
on the wavelength regime (see, e.g., de Jong [19961 , in the sense 
that r ( | is smaller in redder filters. The effect is rather small in 
the F814W filter for the redshifts covered by our sample and at 
maximum corresponds to an overestimate of r& by 11% for the 
highest-redshift galaxy in our sample, at z = 0.97. We corrected 
(i.e. reduced) all measured disk scale lengths for this effect — 
depending on a given galaxy’s redshift — to make them directly 
comparable. 


4.3. Kinematics 

Rotation curves (rotation velocity as a function of radius) were 
extracted from the two-dimensional spectra by fitting Gaussian 
profiles to the emission lines stepwise along the spatial axis. We 
used a boxcar of three pixels, averaging over a given spectral row 
and its two adjacent rows. This approach increases the S /N with¬ 
out a loss in spatial resolution — the boxcar size corresponds 
to 0.6arcsec in the FDF spectra and 0.75arcsec in the WHDF 
spectra; both values are below the average seeing during spec¬ 
troscopy. All detected emission lines were used, and the rotation 
curve with the best S /N and largest spatial extent was used as 
reference in the further analysis. The typical error on the rota¬ 
tion velocity at a given galactocentric radius is 10-20km/s. Ap¬ 
prox. half of the reference rotation curves stem from the [OII] 
doublet, the other half is based on either the [O III], H/i or Her 
line. The vast majority of the rotation curves extracted from dif¬ 
ferent emission lines agreed within the errors of the Gaussian 
fits, in the sense that the rotation velocities at given radius agreed 
within their errors. 

The derivation of the maximum rotation velocity L max of dis¬ 
tant galaxies is a challenging task. At z ~ 0.5, the half-light ra¬ 
dius of a Milky Way-type galaxy is similar to the slit width and 
the typical seeing FWHM. This leads to strong blurring effects 
in the observed rotation curves. It is mandatory to take these ef¬ 
fects into account to avoid underestimates of L max . To tackle this 
problem, we introduced a methodology that simulates all steps of 
the observation process, from the intrinsic 2-D rotation velocity 
field to the extracted 1-D rotation curve. 

The intrinsic rotation velocity field — unaffected by any ge¬ 
ometrical, atmospherical or instrumental effects — is generated 
assuming a linear rise of the rotation velocity V rot (r ) at radii r < 
r t , where r, is the co-called turnover radius, and a convergence 
of Vj-otM to a constant value V max at r > r, (Courteau 1997 ). 
We adopt a turnover radius equal to the scale length r,j <gas of the 
emitting ionized gas disk; rj gas , in turn, is computed from the 
stellar disk scale length r& following Ryder & Dopita ( 119941 ). 
The intrinsic velocity field is then transformed into a simulated 
rotation curve including the following effects: 

1. disk inclination angle ;; 

2. misalignment angle 5 between the slit orientation and the ap¬ 
parent major axis; 

3. seeing during spectroscopy; 

4. luminosity profile weighting perpendicular to the slit direc¬ 
tion; 
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5. blurring effect due to the slit width in direction of disper¬ 
sion — the optical equivalent to “beam smearing" in radio 
observations. 


The simulated rotation curve is then fitted to the observed one to 
infer the intrinsic maximum rotation velocity V max . 

In this paper, we use V max as the only free parameter in the ro¬ 
tation curve fitting process. The other parameters are held fixed, 
based on the observed disk inclination, position angle, scale 
length, etc. For testing purposes, we utilized r t as a second fit¬ 
ting parameter and found that the results on F max agreed with the 
single-parameter fits within the errors for ~ 90% of the objects. 
This is probably due to the strong blurring effects which limit 
or even erase information on the rotation curve shape at small 
galactocentric radii. 

We also investigated whether our results are depending 
on the intrinsic topology of Vj- ot (r). When we, instead of the 
Courteau et al. ( 119971 ) parametrization, use the more complex 
Universal Rotation Curve (URC) shape introduced by Persic et 
al. d 1996k the recomputed U max values agree with the Courteau- 
based ones for 96% of our sample. In brief, the URC was inferred 
from > 1000 observed rotation curves of local spiral galaxies; it 
comprises a mass-dependent velocity gradient. Very low-mass 
spirals show an increasing rotation velocity even at large radii, 
whereas the rotation curves of very high-mass spirals are mod¬ 
erately declining in that regime. The fact that our results on V max 
are not sensitive to the choice of the intrinsic rotation curve shape 
probably has two reasons. Firstly, our sample mainly covers in¬ 
termediate masses, where the URC does not introduce a veloc¬ 
ity gradient at large galactocentric radii. Secondly, the radial ex¬ 
tent of the observed rotation curves — typically four disk scale 
lengths — might be too small to robustly detect a potential ve¬ 
locity gradient in the outer disk. 

To ensure a robust analysis of the scaling relations, we visu¬ 
ally inspected all rotation curves to identify those which a) show 
kinematic perturbations or b) have an insufficient spatial extent 
and do not probe the regime where V mi (r) converges to V max ; in 
extreme cases these curves show solid-body rotation. During this 
visual inspection, we rejected 101 objects from the data set of 
238 galaxies with determined redshifts, corresponding to a frac¬ 
tion of 42 %. Another 13 galaxies were not considered for further 
analysis because the rotation curve fitting yielded a very large 
error on V max , with a relative error cr vre i = cr vmax /V max > 0.5 
(errors on V max stem from fitting the synthetic rotation curves 
to the observed rotation curves via x 2 minimization). The ex¬ 
cluded 13 galaxies are relatively faint (median total apparent R- 
band brightness ( R ) = 22.42, compared to (R) = 21.81 for the 
rest of the sample with derived V max ) and their rotation curves 
show slight asymmetries; otherwise these galaxies have proper¬ 
ties similar to the remaining sample. 

The 124 galaxies with a robust U max constitute our kinematic 
sample. It covers a range 25 km/s < V max < 450 km/s with a me¬ 
dian value of <V max ) = 145 km/s. The median of the relative error 
of the maximum rotation velocities is (<x vre i) = 0.19. This error 
in general is larger towards lower values of V max , i.e. towards 
low-mass galaxies, and smaller towards higher values of Vj nax , 
i.e. towards high-mass galaxies. Table [I] lists the main parame¬ 
ters of ten galaxies as an example (the full table comprising 124 
galaxies is electronically available): object ID, redshift, max¬ 
imum rotation velocity, /(-band absolute magnitude, and disk 
scale length. 


Table 1. Main galaxy parameters 


ID 

z 

^max 

[km/s] 

M B a 

[mag] 

r/ 

[kpc] 

698 

0.5663 

155 ± 46 

-20.43 ± 0.09 

2.90 + 0.58 

745 

0.6986 

290 ± 34 

-21.77 + 0.07 

2.83 + 0.57 

759 

0.0718 

65 ± 15 

-17.14 + 0.09 

1.06 + 0.21 

762 

0.4343 

280 ± 84 

-20.70 ± 0.09 

1.63 + 0.33 

814 

0.6491 

145 ± 55 

-21.09 + 0.09 

4.10 + 0.83 

832 

0.5477 

260 ± 104 

-19.55 ± 0.09 

2.29 ± 0.46 

868 

0.4573 

135 ± 18 

-20.71 ±0.09 

3.18 + 0.64 

870 

0.2775 

90 + 4 

-20.23 ± 0.07 

2.50 ± 0.50 

876 

0.8324 

165 ± 18 

-21.54 + 0.09 

4.68 ± 0.95 

878 

0.2128 

140 ± 50 

-19.09 ± 0.09 

3.65 ± 0.73 


Notes. Ten objects shown as examples. The full table with 124 galaxies 
is electronically available. <a) Rest-frame B-band absolute magnitude, 
corrected for intrinsic absorption following Tully et al. ( 11998 !). <w Disk 
scale length from F814W HST/ACS imaging, corrected for rest-frame 
wavelength dependence following de Jong ( 1996) . 

4.4. The Tully-Fisher relation 

We show the distribution of our kinematic sample in rotation 
velocity - luminosity space in Fig. [2] Our data representing a 
median redshift of z = 0.45 are shown as solid symbols in com¬ 
parison to the local Tully-Fisher relation as given by Tully et 
al. (119981) : 

M b = -7.79 • log Vmax - 2.91, (2) 

which is shown as a dashed line in Fig.0 the dotted lines indi¬ 
cate the 3 cr scatter of the local data. Local and distant galaxies 
are consistently corrected for intrinsic dust absorption following 
the approach of Tully et al. (119981 ). The redshift distribution of 
our kinematic sample of 124 galaxies is shown in Fig. Q] It is 
very similar to the full sample of 238 galaxies with determined 
redshifts. The majority of the distant galaxies are located within 
the 3 cr limits of the local TFR (see Fig. 0. Interestingly, almost 
all galaxies with rotation velocities log V max < 2 fall on the high- 
luminosity side of the local relation, and the eight galaxies with 
log V max < 1.8 are even located above the 3 cr limits of the local 
TFR. Three possible explanations have to be considered. 

Firstly, U max could be underestimated for low-mass disk 
galaxies: optical rotation curves of low-mass disks in the local 
universe often have a positive rotation velocity gradient even at 
the largest covered radii (e.g., Persic et al. 11996 ). As outlined 
in Sect. 14.31 we did not find significant changes in F max when 
adopting intrinsic rotation velocity fields with positive gradients 
in low-mass galaxies (using the prescriptions presented in Persic 
et al.). However, the spatial extent of the rotation curves in our 
sample (as well as other samples at similar redshifts) is typically 
three to four times the optical disk scale length, which might be 
insufficient to constrain potential Vj 0 , gradients. 

Secondly, the observations might hint to a mass-dependent 
evolution in luminosity that is larger for low-mass disk galax¬ 
ies than for high-mass ones. This could be interpreted as a 
manifestation of the down-sizing phenomenon (e.g., Kodama et 
al. 120041 . A third interpretation would be that the apparent mass- 
dependency is only mimicked by the impact of the magnitude 
limit of our sample. It has been demonstrated by numerous au¬ 
thors (e.g. Wi 11 ick l994 ) that, towards lower U max , mainly galax¬ 
ies on the high-luminosity side of a parent unbiased Tully-Fisher 
distribution enter an observed sample. For a detailed discussion 
of this magnitude limit effect, please see Appendix lAl 
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Fig. 2. Rest-frame 6-band Tully-Fisher diagram, showing a comparison between our sample of 124 disk galaxies at a median redshift ( z > ~ 0.5 
(solid symbols) and the local Tully-Fisher relation as given by Tully et al. ( 19981 the black dashed line indicates the fit to the local data (not shown 
in this figure); the dotted lines give the local 3 cr scatter). The blue solid line shows the fit to the distant sample (with the slope fixed to the local 
value), which is offset from the local Tully-Fisher relation towards higher luminosities by (A M B ) = -0.47 ±0.16 mag. 


While the first interpretation outlined above is a purely kine- 
matical one, the other two only concern the luminosity, either 
in the rank of a physical effect (down-sizing) or a selection ef¬ 
fect (magnitude bias). Because our targets have been selected 
for observation using a limit in apparent magnitude — which 
is a common and well-motivated approach — the third scenario 
seems much more likely than the second one. We again refer the 
reader to Appendix lAl for a further discussion of this topic. To¬ 
wards higher F max , however, most rotation curves should be flat 
even beyond the radii probed by our data, provided that the in¬ 
ternal mass distribution in intermediate-z and local disk galaxies 
is similar (e.g. Sofue & Rubin 12001 ). Furthermore, the impact 
of the magnitude limit becomes negligible towards higher V max 
(e.g. Giovanelli et al. [1997 ). 

Note that overestimated luminosities at low T max are highly 
unlikely, as the ground-based photometry is very deep and we 
compute Mg from the filter which best matches the 6-band in 
rest-frame, assuring very small total errors on the absolute mag¬ 
nitudes (combined random and systematic errors are cr MB < 
0.12 mag for all galaxies). The offsets of the slow rotators from 
the local TFR can also not be attributed to an evolution in red- 
shift: the median redshift of these galaxies (( z ) = 0.37 for the 
eight TFR outliers at log V max < 1.8) is lower than that of the 
remaining sample ((z) = 0.45). We will revisit the question of 
any redshift dependencies further below. 

Computing the offsets of individual galaxies from the local 
TFR (as given in Eq. [2j via 

AM b = M b +7.79 -log V max + 2.91, (3) 

we find that the distant galaxies are, for a given V max , more 
luminous than their local counterparts, with a median value 
(A Mb) = -0.47 + 0.16mag (shown as a solid line in Fig. [2}. 
The scatter in A Mb, i.e., the scatter of the distant TFR under 


the assumption of the local slope, is <x 0 b s = 1.28 mag, which is 
2.3 x cr 0 bs of the local 6-band TFR, for which Tully et al. give 
cr obs = 0.55 mag. The fixed-slope fit to the distant data shows 
a significant over-luminosity of the distant disk galaxies. How¬ 
ever, the scatter of the distant sample is large, more than twice 
the local value, and the vast majority of the distant galaxies are 
located within the 3 cr limits of the local TFR. This large scatter 
might be in part due to the broad range in redshifts (see Fig. 0. 
as the evolution in luminosity might depend on look-back time 
(other possible sources of increased scatter will be discussed in 
Sect. 0. We will therefore now shift our focus from the global 
evolution in luminosity to the evolution of individual galaxies. 

We show the offsets from the local TFR, computed following 
Eq. 0 as a function of redshift in Fig. 0 The errors on A Mb 
are computed via error propagation from the errors on F max and 
Mb- The TFR offsets A Mb show an evolution towards higher 
luminosities at higher redshifts. This is confirmed by a linear fit 
to the full sample, which yields 

A M b = -(3.82 ± 1.74)-log(l + z) - (0.15 ± 0.32). (4) 

At given V max , disk galaxies at z = 1 are on average more lumi¬ 
nous by A Mb = —1.2 ± 0.5 mag according to this fit. As shown 
in Fig. 0 the fit to our data is in good agreement with Dutton et 
al. (1201 11 1 who used combined /V-body simulations and semi- 
analytical models. With A-body/hydrodynamical simulations, 
Portinari & Sommer-Larsen (120071) find A Mb = -0.85 mag 
at z = 1; in excellent agreement with the Dutton et al. value of 
A Mb = -0.9 mag, and our own result. The non-linear evolution 
of the model prediction (see Fig. 0 can, however, not be 
confirmed with our data. A second-order polynomial fit agrees 
with the linear one to within +0.1 mag throughout the probed 
redshift range. 
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Fig. 3. Offsets A M b of the galaxies in our sample from the local Tully-Fisher relation, displayed as a function of redshift. The galaxies show 
increasing over-luminosities towards longer look-back times, reaching A M B = -(1.2 ± 0.5) mag at z = 1 according to a linear fit, depicted as 
a blue solid line and dashed lines indicating the 1 tr error range. The observed luminosity evolution is in good agreement with predictions from 
numerical simulations by Dutton et al. ( 2011 i solid magenta line) who find A M B = -0.9 at redshift unity. The dotted line illustrates no evolution 
in luminosity. 


4.5. The velocity-size relation 

We will now investigate the rotation velocity-size relation. For 
this purpose, we utilize the data on ~ 1100 local disk galaxies 
by Haynes et al. (11999b| . Since the electronically available disk 
sizes are given as isophotal diameters at 23.5 mag /-band surface 
brightness, we transformed these into disk scale lengths assum¬ 
ing an average central surface brightness of /r/ = 19.4mag, as 
given in Haynes et al. (11999al) . 

With a bisector fit — a combination of two least-square 
fits with the dependent and independent variables interchanged 
(e.g. Isobe ll990 1) — to the local data, we find 

log r d = (1.35 ± 0.04) • log V max - (2.41 ± 0.08). (5) 

Fig EJshows our data compared to the local VSR as defined 
byEq.0 A few galaxies at the lowest V max are scattered towards 
large radii in VSR space. These are galaxies that also fall above 
the 3 cr-limits of the local TFR. Either V max is underestimated for 
these galaxies, or they are more luminous and larger than to be 
expected for their maximum rotation velocities. 

We compute the offsets from the local VSR as 

A log r d = log r d - 1.35 ■ log V max + 2.41, (6) 

which yields a median of (A log r d ) = -0.10 + 0.05 (displayed as 
a solid line in Fig. The scatter in the distant sample is cr 0 b s = 
0.27 dex; 2.1 x larger than the local value of cr obs = 0.13 dex. 

We now again focus on a potential evolution with look-back 
time. Fig.[5]shows the VSR offsets as a function of redshift. Er¬ 
rors on A log r d are propagated from errors on V max and r d . Using 
a linear fit, we find 

A log r d = -(0.54 ± 0.37) • log(l + z) + (0.00 ± 0.07). (7) 


Disk galaxies at z — 1 therefore were smaller than locally 
by a factor of 1.45 on average, for a given V max . This 
reflects ongoing disk growth with cosmic time, as expected for 
hierarchical structure formation (e.g., Mao, Mo & White 199 8). 
As a test, we also applied a second-order polynomial fit to the 
VSR offsets, not finding a significant second-order term and 
large errors. Similar to the TFR offsets, we therefore can not 
reproduce the non-linear evolution with cosmic time predicted 
by the Dutton et al. (1201 11 1 models. 


5. Discussion 

To interpret the observed evolution of our disk galaxy sample 
in Tully-Fisher space, we need to recall that several processes 
might occur between previous cosmic epochs and the local uni¬ 
verse, and it could be a combination of these processes that gov¬ 
erns the evolution in A Mb- For the following discussion, we 
remind the reader that the maximum rotation velocity V max is 
a proxy for the total (virial) mass of a disk galaxy such that 
M v ir °c V^ ax (e.g., van den Bosch 12002 ). Since the virial mass 
is dominated by dark matter, and the optical luminosity is dom¬ 
inated by stellar light, the optical TFR reflects a fundamental 
interplay between dark and baryonic matter. 

The stellar populations of distant galaxies are likely to have 
a younger mean age; this translates into a lower stellar M/L ra¬ 
tio and, in turn, a lower total (baryonic and dark matter) M/L 
ratio than locally. At given V max , distant galaxies thus would be 
expected to be more luminous than local ones, particularly in 
the /1-band considered here which is sensitive to high-mass stars 
with short lifetimes. The gas mass fractions probably are higher 
towards higher redshifts (e.g., Puech et al. 120101 ). corresponding 
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Fig. 4. Velocity-size relation (VSR): disk scale length r d as a function of maximum rotation velocity V m . dx . We show a comparison between our 
(z) ~ 0.5 data (circles) and the local VSR as found with data from Haynes et al. d 1999bl the black dashed line indicates the fit to the local data (not 
shown in this figure) and the dotted lines depict the 3 tr scatter). The blue solid line shows the fit to the distant sample which is offset to smaller 
disk sizes by (Alogr d ) = -0.10 ± 0.05 dex. 


to higher total M/L ratios and TFR offsets A M B > 0. Further¬ 
more, as less gas has been converted into stars, stellar masses 
M, might be lower at given V max , and this, in turn, would lead 
to lower luminosities and, again, positive TFR offsets. The ob¬ 
servational census on the evolution in stellar mass is somewhat 
unclear at z < 1: e.g., Puech et al. (120081 1 find a decrease in stel¬ 
lar mass of AM, = -0.36^ ^ dex up to z = 0.6, while Miller et 
al. (1201 11 1 give a small and statistically insignificant evolution of 
AM, = -0.04 + 0.07 dex up to z = 1. Using simulations, Dut¬ 
ton et al. (1201 11 1 predict a growth in stellar mass by ~ 0.15 dex 
between z = 1 and z = 0 at fixed U max . 

Note that we do not imply that, for an individual galaxy, the 
evolutionary path in TFR space would at all times be purely in 
luminosity between z = 1 and z = 0. E.g., if a disk galaxy would 
undergo a minor merger, the remnant would, after relaxation of 
the kinematical disturbances, most likely be located at a higher 
fmax and a higher luminosity than before the encounter. During 
or shortly after the minor merger, such a galaxy would not enter 
our kinematic sample, even if it was covered by our observations, 
because it would feature a disturbed rotation curve. This is just 
to give an example that we do not claim disk galaxy evolution 
at z < 1 to solely proceed parallel to the luminosity axis in TFR 
space. A minor merger with a low-mass satellite would proba¬ 
bly lead to an increase of the total M/L ratio, as dwarf galax¬ 
ies have higher dark matter mass fractions than galaxies in the 
M* regime (e.g. Moster et al. 12010 1. Lacking data on gas mass 
fractions, we can not infer dark matter mass fractions for our 
sample. Regarding the ratio M,/M v ; r between stellar mass and 
virial mass, Conselice et al. (120051 1 observationally found a mild 
decrease between z « 1 and z ~ 0, corresponding to a slight in¬ 
crease in total M/L ratio. Using semi-analytic models, Mitchel 
et al. ( 2016 ) inferred a basically constant M,/Mh a i 0 at 0 < z < 1 
(cf. their Fig. 2). 


Out of the processes described above, the younger stellar 
populations towards higher redshifts — potentially in combina¬ 
tion with slightly lower dark matter mass fractions — are likely 
dominating in our sample, since we find an increase in luminos¬ 
ity of A Mg = -1.2 ± 0.5 at given V max towards z = 1, hence 
a decrease in total M/L ratio. Our result agrees with previous 
observational findings e.g. by Bamford et al. (120061 ) or Miller et 
al. (1201 lb . and also with predictions from simulations by Porti- 
nari & Sommer-Larsen (120071 ) and Dutton et al. (1201 It . A stellar 
population modeling of 108 galaxies from our sample (Ferreras 
et al. 120141 ) revealed the well-known downsizing effect: the dis¬ 
tant high-mass disk galaxies began forming their stars at higher 
redshifts and on shorter timescales than the low-mass ones. 

The observed evolution in disk scale length (Fig. [5} reflects 
the growth of disks with ongoing cosmic time. Such an evolu¬ 
tion is expected in an LCDM cosmology with hierarchical struc¬ 
ture growth. Based on theoretical considerations, it has been pre¬ 
dicted already by Mao, Mo & White (119981 ). In the cosmology 
adopted here, the computations by these authors correspond to 
a disk size increase by a factor of ~ 1.9 between z = 1 and 
z = 0 at given V max ; larger than what we find, but almost in 
agreement within the errors. The fit to the observed evolution 
given in Eq. |7] corresponds to Alogra = -0.16 dex at z = 1. 
This is a slightly stronger evolution in size than given in the ob¬ 
servational study of Vergani et al. ( 20121 these authors derived 
A log r ( | = -0.12 dex at z = 1.2), and a slightly smaller evo¬ 
lution than predicted with semi-analytical models by Dutton et 
al. (1201 11 ). who find Alogrd = -0.19 dex at z = 1. Candidate 
processes to explain the disk growth towards z = 0 are accre¬ 
tion of cold gas or minor mergers with small satellites. Any ma¬ 
jor mergers in the cosmic past of galaxies in our kinematical 
sample must have occurred several Gyr ago so that the merger 
remnant could regrow a rotationally supported disk (e.g., Gov- 
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Fig. 5. Offsets A log r A of our galaxy sample from the local velocity-size relation (see Fig. [4j as a function of redshift. Negative values in A log r d 
correspond to smaller sizes at given maximum rotation velocity V max . We find successively smaller disks towards higher redshifts. A linear fit to 
our data (displayed as a blue solid line and dashed lines indicating the 1 cr error range) shows that disk galaxies at z = 1 were smaller than their 
local counterparts by a factor of ~ 1.5. For comparison, the solid magenta line shows predictions from simulations by Dutton et al. 1 2011 1. The 
dotted line illustrates no evolution in size. 


ernato et al. 12009 ) — major merger remnants are kinematically 
cold only in special pre-merger configurations (e.g., Springel & 
Hernquist 12005 ). 

Note that the Haynes et al. d 1999bl> sample that we used as a 
local VSR reference comprises disk scale lengths derived in the 
/-band, the response function of which is very similar to that of 
the F814W filter used for our HST imaging. Since we corrected 
all scale lengths to rest-frame in our data set, local reference 
and distant galaxies can be directly compared. If the required 
correction for the wavelength dependence of r d (following de 
Jong [199 6) would be omitted, our sample would yield a weaker 
size evolution, corresponding to A log r (1 = -O.lldex at z = 1 
for a given V max . 

Our analysis so far has shown that the disk galaxy population 
as a whole is evolving in luminosity and size over the redshifts 
covered by our data. However, we do not know yet how this com¬ 
bined evolution proceeds for individual galaxies. Do the galaxies 
with the strongest evolution in size also show the strongest evo¬ 
lution in luminosity, or is the picture more complex? 

To look further into this, we compare the offsets from the 
Tully-Fisher relation A Mg to the offsets from the velocity-size 
relation Alogrd as shown in Fig. [6] We will first consider the 
situation in the present-day universe: the upper left panel shows 
the local sample from Haynes et al. dl999bt . The offsets A Mg 
of these galaxies are computed from the local TFR as derived 
with the same sample, and also the offsets A log r d from the 
local VSR (as given in Eq. 0 are based on the same sample. 
By construction, the median of both parameters is zero. At first 
glance, it might be surprising that A Mg and A log r d are corre¬ 
lated: the dashed line depicts a fit to the local data; the dotted 
lines illustrate 3 x the local scatter, which in terms of A log r d 
is cr = 0.09 dex. This correlation can be understood as a re¬ 


sult of the fact that disk galaxies populate a plane in the three- 
dimensional parameter space span by maximum rotation veloc¬ 
ity, luminosity and size (e.g., Koda et al. 120001 ). Since the TFR is 
a projection of this fundamental plane, any deviation of a galaxy 
in the local universe from the TFR depends partly on its size, or, 
more precisely, on its position in luminosity-size space. In other 
words. Fig. 6 is equivalent to an edge-on view on the fundamen¬ 
tal plane of disk galaxies in the direction of — but not parallel 
to — the V max axis. 

We divide our sample into three redshift bins z < 0.36, 
0.36 < z < 0.59 and z > 0.59, each holding 41-42 galax¬ 
ies. These sub-samples are shown in the other panels of Fig. [6] 
in comparison to the local fit and scatter. For the distant data, 
TFR offsets A Mg and VSR offsets A log r d were computed as 
described in Sects. I4~4l and l4. 51 To interpret these graphs, it has 
to be kept in mind that any evolution in luminosity and/or size 
is imprinted on the correlation between A Mg and A log r d ex¬ 
plained above. 

We find that, towards higher redshifts, the distant galaxies 
are gradually shifting away from the local AMg-A log r& rela¬ 
tion. Using fixed-slope fits to determine the offsets from the 
local relation in terms of A log r d , we infer an evolution of 
-0.08+0.06 dex, -0.19 ± 0.08 dex and -0.29 ± 0.07 dex at red¬ 
shifts z < 0.36, 0.36 < z < 0.59 and z > 0.59, resp. These offsets 
are depicted by solid lines in each of the distant data panels of 
Fig. [6] These deviations from the local relation are larger than 
the evolution in size alone (see Sect. 1431 since here A log r d is 
a combination of the evolution in size and (in projection) lumi¬ 
nosity. The scatter in the distant AMg-A log r d relation is larger 
than locally: our data yield 0.22 dex at z < 0.36, 0.18 dex at 
0.36 < z < 0.59 and 0.20 dex at z > 0.59. Similar to the situation 
in TFR and VSR space, the scatter in the distant data hence is 
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AMg [mag] AMg [mag] 


Fig. 6. Combined look at the offsets A M B from the Tully-Fisher relation and the offsets A log rj from the velocity-size relation. Top left: local 
sample from Haynes et al. ( 1999b; , no error bars shown for clarity), showing a correlation between the two parameters (the dashed and dotted lines 
show the fit and 3 cr scatter). Top right , bottom left and bottom right: distribution of our data in three redshift bins comprising 41-42 galaxies each, 
with a fixed-slope fit indicated by a solid line. The local fit and 3 a scatter are displayed for comparison in each panel. These figures show the 
combined evolution in luminosity and size towards higher higher redshifts. See text for details. 


approx, doubled compared to the local reference which shows 
cr = 0.09 dex. 

The fact that the correlation between A Mb and A log r ( | holds 
up to redshifts z ~ 1 has consequences for the combined evo¬ 
lution in luminosity and size of individual galaxies. This be¬ 
comes particularly clear in the z > 0.59 bin, which represents 
the longest lookback times (5.6 Gyr < fiookback < 7.7 Gyr) and 
is most sensitive to the evolution with cosmic time. The shape 
of the distribution does not appear to be changed with respect to 
the local universe, but merely shifted towards higher luminosi¬ 
ties and smaller sizes. The galaxies which evolved strongest in 
luminosity are not the ones that evolved strongest in size, and 
vice versa. The galaxies with the strongest decrease in size scat¬ 
ter around A Mb ~ 0, and the galaxies with the strongest evolu¬ 
tion in luminosity mostly show A log rj > 0 and hence are larger 
than their local counterparts at the same V max . 

We finally want to address the scatter of the intermediate- 
redshift scaling relations. The observed scatter at z ~ 0.5 (<x 0 b s = 
1.28 mag and cr 0 t, s = 0.27 dex for the TFR and VSR, respec¬ 
tively) is approx, doubled with respect to the local reference data. 
The TFR scatter we find is smaller than in the sample of Weiner 
et al. (20061 <x 0 b S ~ 1.5 mag), similar to Fernandez Lorenzo et 
al. (2 010 cr 0 b s » 1.2 mag), and slightly larger than the value 
given by Bamford et al. ( 20061 cr 0 b S « 1.0 mag); all studies at 
redshifts similar to our data. Part of the observed distant scatter 
stems from the uncertainties of the V max derivation and the ob¬ 
servational limitations, such as beam smearing, limited spatial 


resolution and so forth. We now want to clarify whether the in¬ 
creased scatter is driven by the measurement errors of the galaxy 
parameters, or due to an increase of the intrinsic scatter. The ob¬ 
served scatter cr obs of the TFR comprises contributions from the 
errors cr m b on the absolute magnitudes, errors cr vmax on the max¬ 
imum rotation velocities and the intrinsic scatter cr,,,, such that: 

^obs — °"mb C Cr vmax + °"int’ 

where c is the TFR slope. Our kinematic sample of 124 galax¬ 
ies yields an intrinsic scatter <Tj nt ~ 1.1 mag. In comparison to 
other studies at similar redshifts, our result is only slightly larger 
than the value given by Bamford et al. ( 2006 ; <x int ~ 0.9 mag), 
but much larger than the <x int ss 0.7 mag derived by Miller et 
al. (1201 Ik Tully et al. (1 1998b do not give the intrinsic scat¬ 
ter for their local TFR sample, but the observed fi-band scatter 
cr 0 bs = 0.55 mag suggests that the intrinsic scatter is of the order 
cr int sb 0.3-0.4mag. Our analysis as well as those of Bamford et 
al. and Miller et al. hence show an increased intrinsic TFR scat¬ 
ter at intermediate redshifts. This increase might be due to, e.g., 
a larger contribution of non-circular motions, i.e. kinematically 
“heated" disks (e.g., Forster-Schreiber et al. 2009 ) or more fre¬ 
quent mismatches between photometric and kinematic position 
angle (e.g. Kutdemir et al. 20101 ) than locally. 

In principle, it would also be possible that the stellar pop¬ 
ulations have an effect on the intrinsic TFR scatter, e.g. due to 
a broader distribution in stellar M/L ratios towards higher red- 
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shifts. To investigate this, we computed the stellar fi-band MjL 
ratios from the absorption-corrected rest-frame (B - R) colors 
following Bell & de Jong (1200 ll) for all galaxies in our sam¬ 
ple. We define three redshift bins, using only the 47 galaxies 
with stellar masses M t > 10 10 M o (which we detect up to red- 
shift z « 1), to minimize the impact of the correlation between 
stellar M/L ratio and galaxy stellar mass. The three resulting 
redshift bins with 15-16 galaxies each have median redshifts of 
( z ) = 0.36, (z) = 0.63 and (z) = 0.84. We find an r.m.s. of the 
fi-band stellar mass-to-light ratio of ctm/l = 0.98 [(M/Lb)q\ in 
the lowest redshift bin, (Tm/l = 0.46 [(M/Lb)o\ at (z) = 0.63 and 
ctm/l = 0.88 [( M/Lb)q ] at (z) = 0.84. At least as far as stellar 
mass-to-light ratios are concerned, we thus find no indication 
that the evolution in the intrinsic TFR scatter might be (partly) 
driven by the stellar populations. 

6. Conclusions 

Utilizing the FORS instruments of the ESO Very Large Tele¬ 
scope, we have constructed a sample of 124 disk galaxies up 
to redshift z ~ 1 with determined maximum rotation velocity 
V ma x- Structural parameters such as disk inclination, scale length 
etc. were derived on HST/ACS images. We analyzed the distant 
Umax - luminosity (Tully-Fisher) and V m - dx - size relations and 
compared them to reference samples in the local universe. Our 
main findings can be summarized as follows: 

1. At given V max , disk galaxies are more luminous (in rest- 
frame B-band) and smaller (in rest-frame /-band) towards 
higher redshifts. By z = 1, we find a brightening of A Mb ~ 
-1.2 mag in absolute B-band magnitude and a decrease in 
size by a factor of ~ 1.5. 

2. The scatter in the Tully-Fisher and velocity-size relations at 
z « 0.5 is increased by a factor of approx, two with respect 
to the local universe. 

3. The observed evolution in luminosity and size over the past 
~ 8 Gyr is in good agreement with predictions from numeri¬ 
cal simulations (e.g., Portinari & Sommer-Larsen [20071 Dut¬ 
ton et al. 201 1). 

4. An analysis of the combined evolution in luminosity and size 
reveals that the galaxies which show the strongest evolution 
towards smaller sizes at z « 1 are not those which feature 
the strongest evolution in luminosity, and vice versa. The 
galaxies with the strongest deviations from the local VSR 
towards smaller disks have luminosities compatible with the 
local TFR, while the galaxies with the strongest evolution in 
luminosity are slightly larger than their local counterparts at 
similar V max . 

In the next paper of this series, we will conduct a comparison 
between the kinematics of distant disk galaxies observed with 2- 
D (slit) and 3-D (integral field unit) spectroscopy, in particular 
with respect to scaling relations like the TFR and VSR (Bohm 
et al., in prep.). Another paper will focus on the evolution of the 
correlation between the maximum rotation velocity of the disk 
and the stellar velocity dispersion in the bulge (Bohm et al., in 
prep.). 
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Appendix A: The distant Tully-Fisher relation slope 

We implicitly assumed in our analysis in Sect. l4.4l that the TFR 
slope remains constant over the redshift range 0 < z < 1. We 
want to justify this in the following. To derive the distant TFR 
slope, we rely on the so-called inverse fit which is of the form 

log Vmax = f(M B ) = a-M b+ b. (A. 1) 

This fitting method is more robust against selection effects aris¬ 
ing from a magnitude limit than the classical forward fit 

Mb = /(Vmax) -C- log ^max "f" d, (A.2) 

as has been demonstrated e.g. by Willick et al. dl995t . Even 
though some statements given in an earlier work of this author 
(Willick Fl994 1t imply that the inverse fitting method is prone to a 
magnitude bias of similar strength as the forward fitting method, 
the results of Willick et al. ( 119951 ) clearly show that the impact 
of the “magnitude bias" arising from sample incompleteness to¬ 
wards lower luminosities is much weaker when an inverse fit is 
utilized. In fact, the strength of the magnitude bias is reduced 
by a factor of six when using the inverse fit, “reducing the bias 
from a significant concern to a marginal effect" (Tully & Cour- 
tois l2012 ). 

Using the parametrization in Ea. IA. II we find 

log V max = (-0.131 + 0.012) • M b - (0.494 ± 0.244) (A.3) 

for our full data set. This corresponds to a distant TFR slope of 
c = -7.62^®g in the form of Eg. IA.2l and agrees well with the 
slope of c — -7.79 for the local sample of Tully et al. ( 119981 ), 
who also used the inverse method. We hence find that the slope 
of the z ~ 0.5 TFR is compatible with the local one. Most other 
observational studies also derived (or assumed) a TFR slope in¬ 
dependent of look-back time (e.g., Bamford et al. 120061 Miller 
et al. 120111 . note that these authors also relied on an inverse fit 
for their analysis). Only Weiner et al. (120061) reported on an in¬ 
creased TFR slope in the range 0.4 < z < 1; however, their 
results are not directly comparable to ours since the approach by 
Weiner et al. lacks corrections for disk inclination and, in turn, 
for the inclination-dependent optical beam smearing effect (see 
Sect.rO. 

In the analysis of an earlier stage of our kinematic survey, 
we found a shallower TFR slope at z » 0.5 using a bisector 
fit (Bohm et al. 2004 ) and showed that an apparent slope evo¬ 
lution in a magnitude-limited survey could be mimicked by a 
strong increase of the TFR scatter with look-back time (Bohm & 
Ziegler[2007]). Applying a bisector fit to our current sample, we 
find a slope of c — -5.02 + 0.47. A forward fit (Ea. IA.21 ) would 
yield an even shallower slope of c — —3.71 ± 0.35. However, 
both these fitting methods are sensitive to the impact of the mag¬ 
nitude bias. To demonstrate this, we use the methodology intro¬ 
duced by Giovanelli et al. (119971 to perform a correction of the 
magnitude bias (as carried out also in Bohm & Ziegler [2007 ). In 
this approach, the observed luminosity distribution is compared 
to a Schechter luminosity function to infer the sample complete¬ 
ness at a given magnitude and, in turn, a given V max . The key 
factor governing the impact of the magnitude bias is the TFR 
scatter, for which we use the observed scatter of our whole sam¬ 
ple, i.e., cr 0 b s = 1.28 mag at z ~ 0.5. After the de-biasing pro¬ 
cedure, the absolute magnitudes are less bright; in particular, for 
galaxies at low V max . For the corrected sample, we find slopes 
of c = -5.82 + 0.35 (forward fit), c = -6.86 + 0.41 (bisector 
fit) and c - -8.33 (inverse fit). These numbers, when com¬ 
pared to the fits of the uncorrected sample, clearly demonstrate 


that the inverse TFR fit is the least sensitive to the influence of 
sample incompleteness. In fact, the inverse fit slopes of corrected 
and uncorrected sample agree within the errors. If we would ap¬ 
ply the de-biasing procedure also to the local sample of Tully et 
al. (which most likely would lead only to very small changes of 
the absolute magnitudes and, in turn, only a very small steepen¬ 
ing of the local TFR slope; however, we could not do this ex¬ 
ercise since the data are not electronically available), local and 
distant inverse-fit slopes would most probably still be in agree¬ 
ment within the fit errors. 

Note that we find large differences between a forward, bisec¬ 
tor and inverse TFR fit not only for our sample but also other 
studies at similar redshifts. Using the sample from Bamford et 
al. (120061) . which comprises 89 field disk galaxies at 0.06 < z < 
1.0, we infer the following TFR slopes: c = -4.59 ± 0.44 (for¬ 
ward), c = -5.89 + 0.56 (bisector) and c = -8.18 ^°^ (inverse). 
Alternatively using the sample from Miller et al. ( 2011 ), which 
contains 129 disk galaxies at 0.20 < z < 1.31, we find these 
slopes: c = -3.99 + 0.45 (forward), c = -5.79 + 0.65 (bisec¬ 
tor) and c = -7.90^074 (inverse). This demonstrates that strong 
differences between the three fitting methods are common for 
distant samples, motivating the use of the method which by far 
is least sensitive to the magnitude bias, i.e., the inverse TFR fit. 
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